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Abstract 

During moonlit nights, observations with ground-based Cherenkov telescopes at very high energies 
(VHE, E > 100 GeV) are constrained since the photomultiplier tubes (PMTs) in the telescope camera are 
extremely sensitive to the background moonlight. Observations with the VERITAS telescopes in the standard 
configuration are performed only with a moon illumination less than 35% of full moon. Since 2012, the 
VERITAS collaboration has implemented a new observing mode under bright moonlight, by either reducing 
the voltage applied to the PMTs (reduced-high-voltage configuration, RHV), or by utilizing UV-transparent 
filters. While these operating modes result in lower sensitivity and increased energy thresholds, the extension 
of the available observing time is useful for monitoring variable sources such as blazars and sources requiring 
spectral measurements at the highest energies. In this paper we report the detection of 7 -ray flaring activity 
from the BE Lac object lES 1727+502 during RHV observations. This detection represents the first evidence 
of VHE variability from this blazar. The integral flux is (1.1 ±0.2) x 10“"cm“^s“' above 250 GeV, which 
is about five times higher than the low-flux state. The detection triggered additional VERITAS observations 
during standard dark-time. Multiwavelength observations with the FLWO 48" telescope, and the Swift and 
Fermi satellites are presented and used to produce the first spectral energy distribution (SED) of this object 
during 7 -ray flaring activity. The SED is then fitted with a standard synchrotron-self-Compton model, placing 
constraints on the properties of the emitting region and of the acceleration mechanism at the origin of the 
relativistic particle population in the jet. 

heat/ingi.'BL Lacertae objects: individual: lES 1727H-502 - galaxies: active - gamma rays: galaxies 
- radiation mechanisms: non-thermal 
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1 . INTRODUCTION 

Blazars are a class of radio-loud active galactic nuclei 
(AON) characterized by a broadband nonthermal continuum 
from radio to 7 -rays, extreme variability and a high degree of 
pol arization. In the framewo rk of the unified AON model (see 
e.g. lUrrv & Padovanilll995h they are considered AON whose 
relativistic jet is aligned with the line of sight. The blazar 
spectral energy distribution (SED) is thus dominated by the 
emission from the jet, enhanced by relativistic effects. In the 
7 -ray sky, blazars are the dominant AON class, representing 
97% of the Fenni-LAJ' AON (between 100 MeV and 100 
GeV, lAckermann et al.l 1201 Ih . and > 90% of extragalactic 
sources detected at very high energies (VHE; £ > 100 GeV) 
by ground-based imaging atmospheric-Cherenkov telescopes 
(lACTsfl 


are indeed HBLs (see e.g. ISenturk et al.ll2013l for a recent 
review). 

The study of VHE blazars is complicated by two observa¬ 
tional characteristics of blazars themselves; their broadband 
emission and their rapid variability. A VHE detection 
on its own, sampling only a small part of the nonthermal 
continuum, does not provide sufficient information about 
the underlying blazar physics. Rather, strictly simultaneous 
multi wave length (MWL) campaigns are required to constrain 
the blazar SED. The launch of the Fermi satellite and its 
monitoring capabilities have had a large impact on blazar 
physics, providing for the first time long-term light curves 
of hundreds of blazars in 7 -rays, and assuring simultaneous 
observations of VHE blazars in the MeV/GeV energy band. 


The blazar class is composed of two subclasses, flat- 
spectrum radio quasars (ESRQ) and BE Lacertae objects, 
depending on the presence (in the former) or absence (in 
the latter) of emission lines in their optical/UV spectrum 
(the threshold between the two subclas ses is an emission 
line e quivalent width of 0.5 nm, see e.g. lAngel & StockmanI 
Il980l) . The two subclasses are also characterized by different 
luminosity and redshift distributions. The ESRQs are on av- 
erage brighter a nd more distant than BE Eac objects (see e.g. 
lPadovanilll992h . This dichotomy in the blazar class reflects 
the dichotomy observed in the radio-galaxy population: ES¬ 
RQs ar e considered the blazar v ersion of the Fanaroff-Riley II 
(FR II, iFanaroff & Rilevl[l974l) radio-galaxies, while BE Eac 
objects are believed to represent the blazar version of FR I. 
The two subclasses share the same SED shape: a broadband 
continuum from radio to gamma-rays, composed of two 
separate bumps, peaking in IR-to-X-rays and MeV-to-TeV, 
respectively. While ESRQs are generally characterized by 
a lower frequency of the first peak, BE Eac objects show a 
vari ety of peak frequenc i es, an d are further classified (see 
e.g. iPadovani & Giommil Il995h into low-frequency-peaked 
BE Eacs (EBE, with Vp^ak in infrared) and high-frequency- 
peaked BE Eacs (HBE, with Vp^ak in UV and beyond). 
When Vpeak is located in the optical/UV the object is often 
classified as an intermediate-frequency-peaked BE Eac (IBE). 

The position of the second peak is related to the position 
of the first one, as shown by observations in the MeV-GeV 
energy band with Fenni-LP^. While ESRQs and EBEs 
present a peak at MeV/GeV, HBEs show a 7 -ray component 
peaking at hi gher energies, ofte n a bove the Fer m i-LP^ 
energy band (lAbdo et'5] 1201 Obh . iFossati et al.l (119981) 
proposed the existence of an anti-correlation between the 
blazar luminosity and Vpeak (the so-called blazar sequence), 
although there is not a gener al consensus on this point. More 
recently iMever et~^ (1201 11) extended this sequence into a 
more general "blazar envelope". The subclass of HBEs, even 
though the least luminous among the other blazar subclasses, 
is the brightest one at VHE, and the majority of VHE blazars 


Observations by lACTs are limited by the high sensi¬ 
tivity of the camera photomultiplier tubes (PMTs), which 
in their standard configuration experience higher noise, 
higher current level, and accelerated aging when operated 
in bright moonlight conditions (moon illumination > 35% 
of full moon). This constraint particularly affects blazar 
studies, for example, limiting the organization of MWE 
campaigns, or prohibiting follow-up observations of a 
flare. Several experiments have successfully performed 
observation s of Cherenkov light under bright moonlight 
conditions (IWeekes et alJ 1 1986 : lOiantell et ^ 11995 1 119971: 

PomarMe et al.l 120011 iKranich et al.l l2001t iTonello et alJ 
20031) . Of the curre nt generation of lACTs, only MAGIC 
and VERITAS (see [Albert et^H] 120071 iRico et all 120081 : 
iBritzger et al.l l2009h observe under moderate moonlight 
(< 35% moon illumination) which significantly improves 
the duty cycle of the observatories. FACT, whose camera is 
composed of solid-state Geiger-mode avalanche photodiodes 
also called silicon-p hotomultipliers), instead of PMTs 
Anderhub et al.l [2013I) . is also capable of observing under 
moonlight. 

In 2012, the VERITAS collaboration began a new program 
of observing under bright moonlight, by applying reduced 
high voltage (RHV) to lower the PMT gain, or by utilizing 
UV-transparent filters. The Schott Glass UG-11 filters’ 
bandpass is from 275 nm to 375 nm, reducing moonlight by 
a factor of ten and Cherenkov radiation by a factor of three. 
The details of the ob serving strategy have been presented a t 
several conferences (iDumm et al.ll2013t IStaszak et al.ll2014h 
and will be discussed in an upcoming VERITAS publication 
(2015, in preparation). In this paper, we concentrate on the 
capabilities of the VERITAS telescope array in the RHV 
configuration. We also present the detection, during May 
2013, of VHE emission from the blazar lES 1727+502 at a 
flux of roug hly five times the ar chival VHE flux measured 
by MAGIC (lAleksic et al.l I20i4l) . The high-flux state was 
initially detected during bright moonlight observations, 
which represents an innovation for VERITAS. 


Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, IL 
60439, USA 

Send off-print requests to: 

Matteo Cerruti (matteo.cerruti@cfa.harvard.edu) 
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* see http: / /tevcat. uchicago . edu for a regularly updated list of 
known TeV sources. 


The b lazar lES 1727+5 02 (from the Einstein Slew Survey 
Catalog. lElvis et al.ll992^■ also k nown asIZw 187 or OT 546, 
is a nearby (z = 0. 055, Okg[1978|) HBE, discovered as a 7 -ray 
source by F ermi (lA bdoet a k 2010al) . and as a VHE source 
by MAGIC (lAleksic et al.l 120141) . Radio observations show 
a compact core-jet rnorphol o gy ( Eaurent-Muehleisen et alj 
1 19931: iKollgaard et alJ 119961 : IPushkarey & Koyalevn20ldk 


typical of blazars, with an apparent jet opening angle of 11 ° 
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jLinford et alJ |2012|). Past optical monitoring shows only 

moderate variability (iPica et al.lll988tlFiorucci & Tostill 19961) 

and t he presence of a weak 300 nm bump (iBregman et al.l 
[1912), quite unusual in BL Lac objects. Prior to being 
detected by MAGIC as a VHE source, lES 1727+502 was 
observed by both HEGRA (lAharonian et al. 120041) and the 
WHIPPLE 10-m telescope (iHoran et al.l2004 ), with no detec- 
tion0 The MAGIC collaboration reported an integral flux of 
(5.5 ± 1.4) X 10“'^ cm“^ s“' above 150 GeV (corresponding 
to (2.1 ± 0.4)% of the Crab Nebula fl u}(Pb and a spectral index 
of 2.7 ± 0.5 tsee lAleksic et al.ll2014l) . 

This paper is organized as follows: in Section |2] we 
describe the capabilities of the VERITAS instrument in the 
RHV configuration; in Section |3] we describe the VERITAS 
observations of lES 1727+502, in both RHV and standard 
configurations; the details of the MWL campaign triggered 
by the detection of the high-flux state are reported in Section 
m in Section 13 we model the SED using a standard one-zone 
synchrotron-self-Compton model, and in Section | 6 ] we 
discuss the results obtained from the MWL campaign. 


2 . y£R/7]45 INSTRUMENT 

VERITAS is an array of four imaging atmospheric 
Cherenkov telescopes located at the Whipple Observatory 
in southern Arizona (31°40’N, 110°57’W) at an altitude of 
1.3 km above sea level. Each telescope is of Davies-Cotton 
design (iDavies & Cottonll 19571) with a 12-m diameter reflec¬ 
tor, and the array is arranged in a diamond configuration with 
^ 100 m to a side. Each reflector comprises 345 identical 
hexagonal mirror facets and has a collection area of 110 m^. 

Each VERITAS telescope is instrumented with a camera 
made up of 499 PMTs each with a 0.15° field of view (EoV), 
for a total EoV diameter of 3.5°. The PMT signals are digi¬ 
tized by 500 megasamples-per-second flash analog-to-digital 
convert ers (EADCs). VERITAS em ploys a 3-level trigger 
system (I Weinsteinll2008l ; IZitzeil20 1 31) and has an array trigger 
rate of ^ 450 Hz. 

VERITAS has an energy resolution of 15%, and a single¬ 
event angular resolution of 0.1° at 1 TeV. During standard 
observations (i.e. when the moon is < 35% illuminated) a 
source with an integrated flux of 1% of the Crab Nebula 
flux can be detected at the 5a level in ^ 25 hours, and a 5% 
Crab source in less than 2 hours. More information on the 
VE RITAS array can b e found in [Holder et al.l (l2006l 1201 Ih 
and iKieda et al.l (120131) . 

In the RHV observation mode, the PMT voltages are 
reduced to 81% of their standard values during dark-sky 
observations, allowing VERITAS to operate when the moon 
is 35-65% illuminated. This reduces the absolute gain 
of each PMT by a factor of ~ 3.2. Observing in RHV 
mode increases the yearly exposure of VERITAS by ^ 13% 
above 250 GeV. Eor RHV data, systematic errors induced 
by uncertainty in the atmosphere, telescope optical point 
spread function, and mirror reflectivity are the same as 

^ A flux upper limit of 8.6% (9%) Crab Nebula units above 300 (940) GeV 
was measured by Whipple (Hegra). 

^ Based on the MAGIC measurement of the Crab Nebula lAlbert et alJ 
12003) 


those for standard data. In terms of sensitivity, a standard 
analysis of VERITAS observations of the Crab Nebula taken 
under dark skies at small zenith angles yields a sensitiv¬ 
ity of ^ 35tT/v^ and an analysis energy threshold of 
^ 170 GeV. In comparison. Crab data taken in RHV mode 
under moonlight also have a sensitivity ^ 35cr/v^ albeit 
with a higher energy threshold (^ 200 GeV). These values 
assume a Crab-like spectral index. Note that VERITAS 
event-selection cuts are optimized for sensitivity and not 
energy threshold. The cuts used in this RHV analysis would 
allow for a 5% Crab source to be detected in less than 2 hours. 


3. VERITAS OBSERVATIONS OE lES 1727+502 

The VERITAS observations that allowed for the detection 
of lES 1727+502 were taken between May 1, 2013 (MJD 
56413) and May 7, 2013 (MJD 56419). There were addi¬ 
tional data taken on May 18, 2013 (MJD 56430) that did not 
result in a detection. After quality selection, approximately 
6 hours of data remain. Of these, 3 hours were taken in 
RHV mode on the first two nights of the exposure. All 
observations were made in “wobble” mode, wherein the 
telescopes were pointed 0.5° away from the target to allow 
simultan eous measuremen ts of the target and background 
regions (iFomin et alJll994l) . and all data in the 2013 dataset 
were taken at elevations between 65° and 70°. 

Two sets of gamma-hadron separation cuts, each optimized 
a priori on data taken on the Crab Nebula where the gamma- 
ray excess has been scaled down to 5%, were used in the 
analysis of the 2013 dataset. For the standard-voltage subset, 
only events with images in three or more telescopes were 
used, whereas in the RHV data subset two-telescope events 
were used and a smaller cut was made on the brightness of 
the shower images to account for the fact that the camera 
gain was reduced. The motivation for using these two sets 
of cuts is to approximately match the energy threshold of the 
analyses, resulting in a threshold of ^ 220 GeV for both the 
RHV and standard voltage datasets. 

The analysis of the complete 2 013 dataset used the 
“reflected-region” background model (lAharonian et al.ll2001l) 
resulting in 159 ON events and 850 OFF events with a back¬ 
ground normalizatio n factor a of 0.0 77, yielding a detection 
significance of 9.3cr (iLi & Malll983L Eq. 17) and an average 
gamma-ray-like event rate of 0.25 ± 0.04 counts/minute with 
a background event rate of 0.18 ± 0.02 counts/minute. 

Prior to the detection of lES 1727+502, VERITAS observed 
the object for 8.6 hours between 2007 and 2009. There was 
no detection in these data; the upper limit (99% confidence 
level) is 2.6 X lO-'^ cm-^s”' (2. 5% Crab) above 350 GeV 
assuming the spectral index from I Aleksic et akl ([20l3)- This 
is consiste nt with the blazar’s low-flux state as detected by 
MAGIC in lAleksic et alJ (120141) . 

The fitted position of the excess events is (a,i 5 )j 2 ooo = 
(17/,28««4.U±7.5^,,t±5.2fy,,-t50° 13'60"±1' 10",,±50",) 
and is within 0.04° (2.4') of the catalog position 
{a,5) = (17^ 2 8"* 18.624",-(-50° 13' 10.416") given in 
iMa et al.l (119981) . The VERITAS catalog name for this source 
is VER J1728-)-502. The excess seen from the target is 
consistent with a point source. 
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TABLE 1 

The data points for the VERITAS observed spectrum of IES 1727+502. 


<E> 

[TeV] 

Elow 

[TeV] 

Ehigh 

[TeV] 

Flux or UL (95% CL) 
[cm“^s“*TeV“*] 

Flux Error 
[cra-2s-‘TeV-‘] 

Excess^ 

Sig. 

0.32 

0.25 

0.40 

3.3 X 10-“ 

9.6 X 10-“ 

25.8 

4.1 

0.50 

0.40 

0.63 

1.1 X 10-“ 

3.1 X 10-“ 

19.2 

4.4 

0.79 

0.63 

1.00 

5.9 X 10-'2 

1.4 X 10-“ 

21.0 

7.0 

1.26 

1.00 

1.59 

1.4 X 10-“ 

5.0 X 10-“ 

8.8 

4.5 

2.00 

1.59 

2.51 

4.1 X 10-“ 

n/a 

0.8 


3.16 

2.51 

3.98 

3.4 X 10-“ 

n/a 

2.0* 


5.01 

3.98 

6.31 

1.6 X 10-“ 

n/a 

1.0* 



^The * indicates that the bin contained no OFF events. 



Fig. 1.— Spectrum for the complete IES 1727+502 dataset. Power-law fits to both the uncoiTected (black, squares) and EBL-c orrected (red, triangles) points 
are provided. The fitted values are given in the text. The blue shaded region represents the result from MAGIC jAleksic et al.l2Ql^ . Note that the first upper limit 
is below the extrapolation of the fit. This may be indicative of a cutoff, however, the number of events in this bin is low (1 ON and 2 OFF events with a = 0.08 at 
2 TeV), thus, it is also possible that this is simply a downward statistical fluctuation. 


The reconstructed VHE spectrum is shown in Fig. □ In 
the same figure we show data points corrected for absorption 
by extragalactic backgro und light (EBL) using the model in 
iFranceschini et al] (120081) (for z=0.055). Both are well fitted 
with a power-law function between 0.25 TeV and 1.6 TeV; 
the x^/NDF for the observed and EBL-corrected spectra are 
1.66/2 {P ^ 44%) and 1.69 jl (P ^ 43%), respectively. The 
observed spectrum is given by 

/ F \ -2.1±0.3 

— = (7.8±l.l)x 10“'U -) cm“2s“'TeV‘, 

dE V620GeVy 

and the EBL-corrected (intrinsic) spectrum is given by 

- = (Ll±0.2)xI0-"(j^) c„,-V'TcV-t 

All errors in the aforementioned fits are statistical. The 
cumulative systematic errors on the flux normalizations and 
spectral indices are conservatively estimated to be 30% and 
±0.3, respectively. This is 50% larger than the standard VER¬ 
ITAS systematic uncertainties. The increase does not have 


a significant impact on the following discussion; a detailed 
discussion of systematic uncertainties during moonlight 
observations will be discussed in the VERITAS moonlight 
paper (2015, in preparation). The numerical values for the 
spectral points are given in Table [1] 

The light curve for the VERITAS observations is shown 
in the top panel of Fig. |2] The peak integrated flux above 
250 GeV is (1.6 ± 0.4) x 10“"cm“^s“', which corresponds 
to 9.5% of the Crab Nebula fluxQ The upper limit is at the 
95% confidence level and represents 9.4 x 10“'^cm“^s“^ 
(5.6% Crab). The light curve between May 01 and May 07 
(MJD 56413-56419) can be fitted with a constant resulting 
in a flux of (1.1 ±0.2) x 10“"cm“^s“* (6.3% Crab) and a 
X^/NDF = 4.73/3 (P - 19%). 

The integral flux above 250 GeV corresponds to about five 
times the flux measured by MAGIC, and represents the first 
evidence of VHE variability in IES 1727+502. The VERITAS 
observations between May 01 and May 07 are consistent with 

^ Based on the power law in IHillas et alj (19981) extrapolated to 250 GeV. 
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a constant flux. However, when including the measurement 
on May 18 in the fit, the resulting -x^/NDF is 12.0/4 (P ^ 
1.8%). Thus, a constant flux is excluded at the 2Aa level, 
indicating that the flare may have ended at some point after 
the last VERITAS detection on May 07. 

4. MULTIWAVELENGTH OBSERVATIONS OE 
lES 1727+502 

4.1. Fermi-LAT 

The blazar lES 172 7+502 was first d etected as a 7 -ray 
source by Fermi-LNT (lAbdo et al.l 120091) . and included in 
both the first and the se cond Fermi catalogs (lAbdo et al.l 
l2010at iNolan et^ 120121) . In the latter, it is named 2FGL 
J1728.2+5015 and is detected with a significance of 
9 .Oct. Its spectrum is parametrized by a power-law func¬ 
tion with index r 2 FGL = 1.83±0.13 and differential flux 
^ 2FCL = (9.5 ± 1.6) X 10“'"^ cm“^ s“' MeV', estimated at 
the decorrelation energy Eq- 2 fcl = 2935 MeV. The source 
has also been included in the catalog of hard Fenni-LALT 
sources (lEHL. lAckermann et al.l2013l) with the name lEHL 
J1728.3H-5014 and an estimated power-law index between 10 
GeV and 500 GeV of Tifhl = 1-67 ±0.34, consistent with 
the 2EGL result. There is no evidence of curvature in its 
Fermi-LAIT spectrum above 10 GeV. 

A new analysis of the Permi-LAT Pass7 (see 
lAckermann et al.l 1201 2^ data has been performed using 
version v9r32p5 of the ScienceTool^ Only photons passing 
the SOURCE class filter and located within a square region of 
side length 40° centered on IFS1727+502 were selected. We 
used photons with energies from 100 MeV to 300 GeV and 
included observations from August 4, 2008 (the beginning 
of the Fermi mission) to August 1, 2013. Eollowing the 
prescriptions of the Permi-LAT Collaboration^ the analysis 
was performed using the binned likelihood method, and the 
P7REP_S0URCE_V15 instrumental response function. Data 
were filtered considering only zenith angles lower than 100 ° 
and rocking angles lower than 52°. The Galactic diffuse 
component and the isotropic background were modeled 
using the templates provided by the Fermi team, namely 
gll_iem_v05.fits and iso_source_v05.txt. The likelihood 
analysis was performed using gtlike and all the 2PGL sources 
present in the region of interest were included in the model, 
using the script make2FGFxml.pyW Eor every source the 
spectral function (log parabola or power law) used in the 
2PGL catalog was adopted, with parameters free to vary. The 
model also included sources outside the region of interest, 
up to a distance of 35°, with spectral parameters frozen to 
the 2PGL values. The normalizations of the models for the 
Galactic diffuse component and the isotropic background 
were left free to vary during the fit. 

Using this dataset, IFS 1727+502 was detected with a test- 
statistic (as defined in lAbdo et al.ll2010ah of 219, correspond¬ 
ing to a significance of ^ 14 ct. Its emission is parametrized 
by a power-law function with index To.i-joocev = L91 ± 0.08 

^ See http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
software/ 

^ See http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
LAT_caveats.html 

See http: //fermi .gsfc.nasa.gov/ssc/data/analysis/ 
user/ 


and differential flux $o.i- 300 GcV = ( 2.1 ± 0 . 2 ) x 10 “'^ 
cm“^ s“' MeV“', estimated at the decorrelation energy 
Eo; 0 . i- 300 GcV = 2136 MeV, consistent with the 2PGL result. 
The presence of curvature in the spectrum was tested by 
replacing the power-law model with a log-parabolic one; 
the value of the curvature parameter /3 is consistent with 
zero, supporting the power-law assumption. In the residual 
map there is no evidence of additional sources besides the 
ones included in the 2EGL catalog. Spectral points were 
computed for five different spectral bins, using the python 
scripts prepared by the Fermi-LAT Collaboration0 assuming 
a fixed spectral index ro,i- 3 ooGcV = L9L If the significance is 
lower than 5 ct, 95% upper limits are provided. The bow-tie of 
the best-fit Fermi measurement, as well as the spectral data, 
are plotted in Eig. [3 

The 2EGL cata log reports a variability index of 18 (see 
INolan et'n]l2012l for details), indicating a lack of variability. 
A test for variability in the MeV-GeV band has been per¬ 
formed using aperture photometry: only photons detected in 
a l°-radius region around the blazar have been considered, 
and count-rate light curves have been produced using gtbin 
with time bins of one day, one week, one month and three 
months. All four light curves are consistent with a constant 
flux. 

A second, unbinned, analysis was performed in order to 
produce a measurement in the GeV energy band spanning 
the same interval as VERITAS observations. Only Eermi-LAT 
observations taken between May 1 and May 7, 2013 (MJD 
56413 and 56419), inclusive, are considered. In this case 
the region of interest is a circle of 20 ° radius centered on 
IFS 1727+502. The model includes all the sources used 
in the binned analysis of five years of Fermi-LPTY data, 
with spectral parameters frozen if the source is outside 
the region of interest. IFS 1727+502 is not detected by 
Fermi-LPTT in this very short period and only upper limits 
on its GeV emission can be computed. They are esti¬ 
mated at the 95% confidence level for three different energy 
bands (0.1-1; 1-10 and 10-300 GeV), and are plotted in Eig. [3 


4.2. Swift-XRT 

The detection of the high-flux state at VHE by VERITAS 
triggered X-ray and UV observations by the Swift satellite 
(iGHirels et al.ll200^ . which observed IFS 1727+502 on May 
7 and May 8 , 2013 (MJD 56419 and 56420) for a total live 
time of 6.8 ks. The X-Ray Telescope (XRT, iBurrows et al.l 
120051) performed observations in windowed-timing mode, 
where only the 200 central columns of the detector are 
read, improving the t ime resolution of the instrument (see 
ICusumano et al.ll2012^ . The data analysis is performed using 
HFASoft (version 6.13). Cleaned event files are produced 
using default screening criteria. Images, light curves, and 
spectra are extracted (using XSelect, version 2.4b) from a box 
with height equal to 10 bins and length equal to 40 bins for 
both the source and the background region. 

The XRT light curve above 0.3 keV (in counts per second), 
corrected for the exposur^B and for background, is shown in 

* See http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
scitools/python_tutorial.html 

^ See http://www.swift.ac.uk/analysis/xrt/Iccorr. 
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Fig. 2.— Light curve of lES 1727+502 during the May 2013 multiwave¬ 
length campaign. The time axis is in Modified Julian Days. For reference, 
MJD 56415 con'esponds to May 3, 2013. Top panel: VERITAS flux light 
curve, measured above 250 GeV and in daily bins. The errors are statisti¬ 
cal only, and the first two data points were taken in RHV mode. The red 
line corresponds to the average VERITAS flux, while the dotted black line 
corresponds to the archival VHE flux reported by MAGIC. Central panel: 
Swift-XRT light curve, expressed in counts per second above 0.3 keV. The 
red line con'esponds to the average X-ray flux during this campaign. Note 
that in this panel the X-axis is zoomed on MJDs 56419 and 56420. Bottom 
panel: 5w/yf-UVOT (on MJD 56419-20) and FLWO 48" (on MJD 56430) 
light curve, expressed as lyF^. Data have been corrected for both Galactic 
absoi'ption and host-galaxy contamination. Red points correspond to i’ and 
r’ filters (the two are almost superposed on MJD 56430), orange to V, blue to 
B, violet to U, magenta to the UVW1,UVW2 and UVM2. 


Fig. m No significant variability is detected within a single 
observation, nor between the two observations: a fit of the 
light curve with a constant function yields a x^/NDF value 
of 714/682 (chance probability of 19%), consistent with the 
flux being steady (within the statistical uncertainties). The 
mean count rate is 0.83 ±0.12 counts per second. The rate 
is low enough to avoid any significant pile-up effect in the 
detector. The data shown in Fig. |2]have been rebinned for 
plotting purposes. 

The spectral analysis is performed using XSpec (version 
12.8.0). Given the lack of variability, and in order to improve 
the statistics, the data are summed (using mathpha) and 
rebinned (using grppha) assuming a minimum of 50 counts 

php 


per bin. The response files provided by the Swift science team 
are used, while the ancillary response files are computed 
using xrtmkarf. Data below 0.3 keV are excluded, and the 
last significant spectral bin extends up to 8 keV. 

The first spectral model tested is a simple absorbed power 
lawPl. The neutral absorption is fixed to the Galactic value 
Nh = 2 .75 X 10^° cm“^, as measured by [Dickev & LockmanI 
The best-fit result is F = 2.22 ± 0.04, with a 
normalization factor C = (4.9 ±0.1) x 10“^ cm“^ s“' keV“', 
estimated at 1 keV. However, the /NDF is 108/86, and 
significant residuals are seen above 5 keV. The fit is sig¬ 
nificantly improved (F-test probability equal to 8 x 10“^) 
if a break is added to the model, considering an absorbed 
broken-power-law function. The best-fit parameter values 
are Fi = 2.0l!;j:|°, Fa = Ebreak = keV, 

and normalization C = (5. 3 ^ 02 ) ^ 10”^ cm“^ s“* keV“'. The 
X^/NDF value is 86/84. The measured (i.e. not corrected 
for Galactic absorption) flux in the energy range 0.3-8 keV is 
Fx = (2.06!“;*]^) X lO”" erg cm-^ s”'. 

The observed spectral break could be intrinsic or due to 
an underestimation of the neutral absorption: in particular, 
additional absorption due to the blazar host ga laxy should be 
studied (for a detailed ana lysis of this effect see lPerlman et al.l 
l2005HFurniss et al.ll2013h . This hypothesis is tested by fitting 
the data with an absorbed power-law function, letting the 
yalue of Nh to yary. The amount of absorption required 
to mimic the break is Nh = (b.Silij'y) x 10^° cm“^, which 
is about 2.4 times the Galactic absorption estimated by 
!Dickey & Lockm^ (Il990h . The yf /NDF yalue of this fit is 
howeyer only 95/85 (chance probability of 0.21; the F-test 
probability between this model and the broken-power-law fit 
is equal to 4 x 10“^), suggesting that the intrinsic origin of the 
break is preferred. The spectrum shown in Fig. |3]is the one 
computed assuming a broken-p ower-law function, correcte d 
for the Galactic absorption from iDickey & LockmanI (1 19901) . 

Swift-XKY also obseryed lES 1727+502 during 201CR 
and th ose obseryations were an alyzed by both iMassaro et al.l 
dMl and lAleksic et al.l (120141) . The count rate during these 
preyious obseryations is lower (around 0.6 counts per second, 
73% of the 2013 rate), suggesting that the VHE high-flux 
state seen by VER ITAS was associated with a higher-flux 
X-ray state. While lAl eksic et aT1 (120141) performed a simple 
power-law fit, IMassaro et al.l (1201 ih confirmed the deyiation 
from a power law, successfully fitting the Swift-XRT data 
with a log-parabolic function. 

4.3. Swift-UVOT 

The UVOT telescope (iRoming et al.l 120051) . on board 
the Swift satellite, obseryed lES 1727+502 at optical and 
ultrayiolet wayelengths, simultaneously with XRT. All mea¬ 
surements were performed using the six ayailable filters: V 

Ttie Galactic absorp tion is computed using the tbnew model, an 
updated version of tbabs IWilms et allI2000h . See http ://pulsar. 
sternwarte.uni-erlangen.de/wilms/research/tbabs/ 

See http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/ 
w3nh/w3nh.pi 

See as well the automatic Swift-XRT preliminary analysis, which in¬ 
cludes the two observations presented here: http: / /www. swift. psu . 
edu/monitoring/ 
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and B in the optical, and U, UVWl, UVM2 and UVW2 in the 
ultraviolet. A circular aperture with radius 5" is used for the 
source, while the background is estimated from a larger re¬ 
gion with radius 15". The flux is estimated using uvotmaghist 
(version 1.1). As in the X-ray observations, no variability was 
detected either within or between the two UVOT observations. 

The Galactic extinction is taken into account assum¬ 
ing Eb-v = 0.037, a value cons istent with the Nh valu e 
used for the X-ray analy sis (iJenkins & Savaga [1973). 
ISchlaflv & Finkbeined (1201 Ih estimated a similar value of 
Eb-v = 0.026. Assuming Eb-v = 0.026 instead of 0.037 
affects the flux estimations by between 3% and 9%: this 
additional systematic uncertainty is included in the UVOT 
error bars plotted in Fig. The correct i on fac tor for each 
Alter is computed followin mRoming et al.l (l2009h . 

In infrared and optical light, the SED of lES 1727+502 is 
dominated by the host galaxy , a giant elliptical well-studied 
in the past (IScama et al.1l200(ilUiTv et aOioOPtlNilsson et al.l 
12001 . For the purpose of subtracting the host-galaxy 
contamination, we made use of the recent results from 
iNilsson et al.l (120071) . who estimate a contaminating flux of 
1.01 ±0.06 mJy in the R band, for an aperture radius of 
5". This contribution is then tra nslated into a correcti on in 
the V, B and U Alters following [Hvvonen et al.l (120071) . An 
uncertainty in the host-galaxy color of 0.1 is considered, 
and included in the Anal error bars of the UVOT spectral 
points. For the remaining UVOT Alters we assume that the 
host-galaxy contribution is negligible compared to the AGN. 


4.4. FLW0 48" 

As part of a long-term optical program of monitoring of 
VHE blazars, lES 1727+502 is regularly observed by the 
automatic FLWO 48" telescope, located near the VERITAS 
siteQ The instrument, a reflector of 1.2-m diameter, can 
perform measurements using the standard Cousins and 
SDSS Alters. The analysis uses an aperture radius of 10" to 
estimate the source magnitude. There are no observations 
performed simultaneously with VERITAS: the observation 
closest to the VERITAS detection of lES 1727+502 was 
performed on May 18, 2013, eleven (ten) days after the last 
detection by VERITAS (Swift), using B, r’, and i’ Alters. 
The data are dereddened using Eb-v = 0.037, as for UVOT, 
and the host-galax y contribution is subtracted following 
INilsson e t al.l (120 071) and t ranslated into the rel evant Alters us¬ 
ing [ Hvvonen et alT(i2007l) (for the B Alter) or iFukugita et al.l 
(II 995 I) (for the remaining SDSS Alters). Similar to UVOT, 
the uncertainty in the dereddening and the host-galaxy color 
is taken into account and included in the error bars. The three 
FLWO 48" measurements, corrected for both absorption and 
host-galaxy contamination, are plotted in Fig. |3] 


5. SED MODELING 

The SED of lES 1727+502 is shown in Fig. |3l Si¬ 
multaneous observations (with Swift-XRT, Swift-UVOT, 
and VERITAS) were carried out only during MJD 56419. 
Given that none of the instruments involved in the campaign 
detected any signiflcant variability, we make the assumption 

http://www.sao.arizona.edu/FLWO/48/48.html 


TABLE 2 

Parameters used for the SSC modeling of IES 1727+502. 


Parameter Value 


5 30 


7mm [103] 

1-5 

y break [10^1 

2.2-3.2 

"/max [to®] 

5.5-7.0 


2.2 


3.2 

K), [10-'* cm-"] 

3.5-29.3 

*Ue [10“^ erg cm"^] 

1.6-12.8 

*Upi, [IQ-® ergcm'3; 

1 8.5-41.5 

Rsrc [10*’ cm] 

4.3-7.4 

*Tyar[days] 

5.8-10.0 

B [mG] 

0.3-0.6 

*ub [IQ-* erg cm"^] 

3.6-16.2 

*Ue/uB [10^] 

1.6-13.8 

*Upi,luB 

1.2-2.8 

*Lje, [10‘'"'erg s'*] 

1.4-5.3 

Note. — Italicized values are fixed in the fitting algorithm. Derived quan¬ 
tities are marked with a star. Two different values of 'fmin are studied, and do 
not affect the other free parameters. The values of 'fmax are simply equal to 

22 X ‘y break- 


that the flux was constant during the observations, and we 
thus consider that the average spectral measurements from 
the different instruments are representative of the blazar 

emission. It is possible that the 

emission was variable, but 


at a level below instrumental sensitivities. Nonsimultaneous 
FLWO 48" measurements also cover the infrared-to-optical 
part of the SED. In this case as well, given the absence of 
variability in the B Alter between Swift-\]WOT and FLWO 
48" observations, we make the assumption that the FLWO 
48" measurements can be considered as representative of the 
emission from IES 1727+502 during the VERITAS observa¬ 
tions. Fenni-LAI upper limits calculated from observations 
performed during the VERITAS campaign (but not strictly 
simultaneously) are also included, as well as the average 
spectrum from the long-term analysis. The bow-tie from the 
MAGIC detection and the archival MWL data from the NED 
(NASA/IPAC Extragalactic Database^ are also included. 

In the context of the synchrotron-self-Compton (SSC) 
model, the two components of the blazar SED are associated 
respectively with synchrotron emission from leptons (e^) and 
inverse-Compton scattering of the particles off synchrotron 
radiation from the same lepton population. The emitting 
region is a spherical blob of plasma (characterized by its 
radius R) in the relativistic jet, moving towards the observer 
with Doppler factor 5, and Ailed with a tangled, homogeneous 
magnetic held B. The particle population is parametrized 
by a broken power-law function, and it carries six free 
parameters; the minimum, maximum, and break Lorentz 
factors of the particles Itmuc, 'Ytreak), the two indices 

(ai and a 2 ), and the normalization factor K. The minimum 
Lorentz factor of the leptons can be Axed at a reasonably 
low value without affecting the modeling. The remaining 
eight free parameters can b e constrained by observations , 
as discussed, for exam ple, in [Bednarek & Protheroel (119971) : 
iTavecchio et alJ (119981) and ICerruti et al.l (l2013al) . 

http: / /ned. ipac .caltech.edu 
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Fig. 3.— Quasi-simultaneous SED of lES 1727+502 (in blue; from low to high energies: FLWO 48", Swift-XKt, Fermi-hAI and VFRITAS). 

Archival data (from the NED, Fermi-hPH, and MAGIC) are included, and plotted in gray. The light-blue and pink curves represent the two sets of SSC models 
computed assuming 7,„,„ = 10^ and 5 X 10^ (see Table|2]and Section|5). The black model represe nts the solution which min imizes the value. The black dotted 
line represents the same model without absorption on the EBL, estimated using the template bv IFranceschini et al.l j2008f) . The luminosity (y-axis at the right) 
is computed assuming z = 0. 055 and the cosmological values (flat universe, with Hq = 67.3 km s“* Mpc“* and Q„, = 0.315) estimated using the measurements 
made by the Planck satellite jPlanck Collaboration et al.l2014l) . 


Given the lack of a simultaneous Fermi-LAT detection, it 
is impossible to constrain the position and the luminosity of 
the inverse-Compton peak during the high state, and a unique 
solution for the SSC model cannot be provided. However, the 
synchrotron component is very well sampled, and can provide 
some constraints on the energy distribution of particles in 
the emitting region. At low energies, the subtraction of 
the host-galaxy contamination reveals the AGN nonthermal 
continuum, which can be described by a power law from 
infrared to UV, with no sign of break. A fit of the FLWO 
48" and UVOT data results in an index ni = 1.64 ±0.09. In 
the case of synchrotron radiation, this index reflects directly 
the index of the underlying e=*= population a\ = 2ni - 1, 
which is thus equal to 2.28 ±0.18. Similarly, the index of 
the X-ray power law below the X-ray break corresponds to 
a 2 = 3.0±0.2. In the following, a\ has been hxed to 2.2 and 
a 2 to 3.2, i.e. a spectral break of 1.0 has been assumed. 

In blazar physics, an additional observational constraint 
comes from the variability timescale r^ar and the causality 
argument: the emitting region size R has to be smaller than 
c7Var^/(l -bz). However, for the case of lES 1727+502, there 
is no evidence of variability on short timescales in any of the 
light curves that can provide an estimation of T,,ar and thus 
constrain the emitting region size. 

As discussed above, the value of 7 m,„ does not affect 
the modeling and it is therefore held fixed. However, 
values lower than 10^ would overestimate archival radio 
measurements (see Figure [3]l, even though, given their non¬ 
simultaneity, these data should not be considered as a strong 


constraint. An additional constraint on is provided by the 
Fermi-LAT nondetection, and will be discussed in the next 
section. In the following we study two different cases, for 
7 mm = 10^ and 5 x 10^. The value of 'y^ax is constrained by 
the break observed in the X-ray spectrum of lES 1727+502. 
In the following, we express it as a function of 'ybreak- 
Imax = 22 jbreak- The numerical factor corresponds to the 
ratio between the X-ray break and the synchrotron peak, 
which is estimated by extrapolating the UV and X-ray spectra. 

Once the indices and maximum energy of the particle 
population have been fixed, the number of free parameters 
in the SSC model is thus five (5, B, R, 'ybreak, and K), and, 
with only four observables (the frequency and flux of the 
synchrotron peak, the VHE spectral index and the VHE flux 
at 620 GeV: I'syn-peak', vFy-syn-peak', ^VERITAS and vF^-yERITAs), 
a unique solution cannot be provided. However, assuming a 
reasonable value of the Doppler factor, it is possible to study 
the parameter space of the remaining free parameters. The 
best-fit solut ion is compu t ed usin g the numerical algorithm 
described in ICerruti et al.l dloTsH)- The parameter space is 
systematically sampled, producing a set of SSC models and 
computing for each of them the expected observable values. 
A fit is then performed in order to express each observable as 
a function of the free parameters, defining a set of equations 
which is solved for the particular set of observables of 
lES 1727+502. The uncertainty in the observables is taken 
into account by iteratively solving the system of equations, 
producing a family of solutions and determining the allowed 
range for each free parameter. 

Eor the particular case of lES 1727+502, the Doppler 
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factor has been fixed to 30 (a value in lin e with standard 
SSC modeling of VHE blazars, see e.g. iTavecchio et alJ 
I 2 OIOI: IZhang et'afl |2014|) . while the other free parame¬ 
ters have been studied in the range: B € [0.1,2] mG; 

e [6 X 10'®, 10^*^] cm; K' € [5 x 10-'2,5 x lO-'"] cm-^; 
'Jbreak ^ [10®, 10®]. K' is defined as the particle den¬ 
sity at 'ybreak- K' = 'The solutions of the SSC 

model are computed iterating on the four observables; 
V„ak e [10'®-*',10'®-«®]; vF,,syn-peak £ [lO-'''''^, lO-'O®^]; 
vF,, VERITAS £ [10-"-4®,10-" '*] and Tveritas £ [1.78,2.52], 
For the two VERITAS observables, the systematic uncertainty 
has been taken into account and summed in quadrature with 
the statistical error. The system of equations is completed by 
an inequality relating the variability timescale to R and 5. We 
considered T^ar =10 days, which corresponds roughly to the 
interval between the last VERITAS detection on MJD 56419 
and the observation during MJD 56430, which indicated 
that the VHE flare may have ended. All the SSC models 
that correctly describe the lES 1727+502 SED are then 
recomputed, and plotted in Fig. |3] The parameter values are 
listed in Table |2] together with derived quantities such as the 
energy densities, the luminosity of the emitting region and 
the minimum variability timescale. For each solution, the 
with respect to the observational data is calculated in order to 
find the solution that minimizes the x^. 


6 . DISCUSSION 

The VERITAS detection of VHE 7 -ray emission from 
lES 1727+502 triggered MWL observations that allowed, 
for the first time, the study of the SED during a high 7 -ray 
flux state, a factor of five higher than the MAGIC detection. 
Even though the lack of simultaneous detection by Fermi- 
LAT does not allow the parameter space to be fully 
constrained, the available data permit some important con¬ 
clusions to be drawn regarding the properties of the particle 
population and the acceleration/cooling mechanisms. 

In the presence of synchrotron radiation, a break in the 
stationary particle popul ation is expected, and it is character¬ 
ized by Aa = 1 (see e.g. llnoue & Takaharall 19961) . Deviations 
from Aa = 1 can, however, occur if the emitting region is 
inhomogeneous, or if the e mission from escaped particles is 
taken into account (see e.g. ISokolov et al.ll20()^ IChen et aP 
l20L5h . Interestingly, the measurements of the synchrotron 
component of lES 1727+502 are consistent with a simple 
break of 1 . 0 , indicating that a single particle population, 
injected with a power-law distribution with index equal to 
2.28 ±0.18, could explain the blazar emission. 

The break energy corresponds to the equality of the 
escape (or adiabatic) timescale and the synchrotron and 
inverse-Compton cooling timescales. The former is energy 
independent (if the escape is advective, and not diffusive), 
and equal Xo R/ /3c, while the cooling timesca l es are both pro¬ 
portional to 1 / 7 . Following ITavecchio et al.l (1 19981 equation 
30), it is possible to relate Xhreak to the model parameters 
and find the only unknown variable, /3 £ [1 /225,1 /70]. This 
means that the stationary particle population is consistent 
with an injected power-law particle population cooled by 
synchrotron and inverse-Compton emission only if the 
advective escape timescale is ~ 1 00/?/c. This v a lue is 
low compared to the ones used by ITavecchio et alJ (Il998l) 


(/3 £ [0.33,1]). For example, assuming /3 = 1, the break 
energy would have been at Xbreak 10 ^, two orders of 
magnitude higher than the constraint provided by quasi- 
simultaneous observations. Alternatively, the hypothesis of 
a power-law population injected and cooled by synchrotron 
and inverse-Compton processes could be too simplistic, and 
additional injection/cooling/escape terms could contribute to 
the final stationary particle population. 

It is important to recall that these considerations are valid 
only for a given value of (5 = 30. Lower values of 5 would 
imply higher values of B and then higher values of /3. For 
example, assuming 5* = 15, the overall normalization is 
reduced by a factor of 2 ^' = 16 and can be compensated by 
assuming a magnetic field four times stronger, B* ^ 1.6 mG. 
In order to maintain the same i^syn-peak (which is oc 
Ibreak has to be reduced by a factor of s/2. These transfor¬ 
mations would affect the inverse-Compton peak as well, and 
it is not guaranteed that these parameters could reproduce 
the VHE spectral index. Regardless, using these new values, 
P* oc Ibreak^*^ — 11/3, One Order of magnitude higher than 
that computed for S = 30. 

It is possible to extract information on the acceleration 
mechanism from the value of the index of the particle 
population. Acceleration by diffusive shocks naturally 
produces power-law particle populations, with an index 
close to 2.0. For ultrar elativistic shocks, several authors (see 
lAchterberg et al.l[ 200 lL and references therein) have shown 
that ai is expected to be around 2.2. The quasi-simultaneous 
multiwavelength (from infrared to UV) observations of 
lES 1727+502 show that ai has to be between 2.10 and 2.46, 
fully consistent with ultrarelativistic shock acceleration. 

The maximum particle energy is another important 
model parameter related to the acceleration. The break 
observed in X-rays at ~ 1.2 keV can indeed be related to 
the high-energy cutoff of the particle population {'fmax), and 
can thus be used to constrain the maximum energy of the 
leptons. In the framework of the standard diffusive shock 
acceleration, the maximum particle energy can be computed 
assuming that the acceleration timescale equals that for 
the radiative losses; particles with energy higher than j,„ax 
lose energy faster than they can be accelerated. Assuming 
Tacc = V mcxIeB = (5.7 x 10®s) "qj/B, where q represents a 
parameter characteristic of the shock acceleration, it is easy 
to estimate j,„ax = 3.7 x 10® {qB„tc)~^'^, which is consistent 
with the constrained values of 7 ,„a^ only for q of the order 
of 10®. Such a high value, which implies rather inefficient 
acceleration, is t ypical of high-frequency -peaked BL Lac 
objects (see e.g. Ilnoue & Takaharal Il996^ . In this case, 
solutions computed for lower values of 6 cannot significantly 
lower the value of q. It is easy to show that, using the same 
values estimated in the discussion on qbreak, q would be 
reduced only by a factor of two. 

The minimum particle energy also carries important 
information on the physics of the emitting region. Both 
archival radio measurements and the nondetection by Fermi- 
LAT suggest that has to be at least of the order of 10®. 
The study of the models for two different values of 7 „„„ 
(equal to 1 and 5 x 10®) shows that solutions with lower 
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values of 7 „„„ have higher emission in the 0.1-1 GeV energy 
band. A 7 mm lower than 10^ would start conflicting with 
the Fermi-hP^ nondetection (under the assumption that the 
VERITAS flux is representative of the average flux throughout 
the whole period used to compute the Fermi-LAil upper 
limit s). Such a high v alue of 7 mm is not a surprise for HBLs 
(see lAliu et aO l2014h . It can be explained, for example, 
by assuming that the particles are already injected with a 
truncated-power-law distribution, and they do not have time 
to cool down completely, or that the parti cles are somehow 
reacc elerated inside the blazar region (iKatarzvhski et alJ 

[200i). 

It is also interesting to note that the one-zone SSC modeling 
suggests that the blazar GeV spectrum during the high-flux 
state was harder compared to the average spectrum measured 
by ferm/-LAT. A harder-when-brighter correlation in the 
Ferm/-LAT band is not common for HBLs (see for example 
lAbdo et aP 1201 lalbh . and is more similar, for example, to 
what is observed during the 7 -ray flaring activity of BL 
Lacertae (lArlen et al.ll2013h . 

One of the open questions in blazar physics is related 
to the energy budget of the emitting region. Is it close to 
equipartition between particle, photon, and magnetic energy 
densities («(., Upi, and ub)1 The solutions estimated here 
are all characterized by a very high value of UejuB, which 
is of the order of 10^“^. On the other hand, the photon 
energy density (synchrotron emission by primary leptons) 
is of the same order as the magnetic one. Equipartition is 
sometimes useful to constrain the parameter space of blazar 
models, especially when dealing with external photon fields, 
which increase the number of free parameters. Recently 
iDermer et al.l(l2Q14ll presented a general approach to compute 
blazar models close to equipartition, showing that, within 
this framework, it is possible to describe the SED of blazars, 
and in particular ESRQs, correctly fitting the spectral break 
observe d by ferm/-LAT in several low-frequency-peaked 
blazars (ICerruti et al.ll201^ . On the other hand, it is not 
clear if a general equipartition approach still holds for HBLs, 
and indeed several modeling attempts (within a standard 
one-zone SSC scenario) on the most well stu died HBLs 
result in equipartition factors far from unity (see lAbdo et al.l 
I 2 OI lallbtlAbramowski et al.ll 2012 h . 

It is important to emphasize that the SED of lES 1727+502 
represents a snapshot of a high-flux state, and that the ratio 
between the flaring and the archival fluxes is higher in the 
inverse-Compton component than in the synchrotron one. It 
is thus possible that the emitting region is only temporarily 
far from equipartition, and that the low state is indeed char¬ 
acterized by Ue/uB closer to one. To test this hypothesis, the 
constraining algorithm has been rerun assuming vF^yERiTAS 
lower by a factor of five (i.e. consistent with the MAGIC 
detection). As expected, the equipartition factor is lowered 
by one order of magnitude, lying between 50 and 1000. 
Again, solutions computed for lower values of 5, implying 
a higher value of B, would result in Ug/uB much closer to 
equipartition. In addition, since the particle index ai is 
softer than 2.0, as constrained by optical/UV observations, 
the value of is strongly dependent on -fmm- Given that we 
cannot really constrain the value of 7 mm, it is possible that the 
equipartition factor is indeed closer to unity. 


7. CONCLUSIONS 

In this paper we presented the results of a multiwavelength 
campaign carried out on the BL Lac object lES 1727+502 
during May 2013, triggered by a VHE high-flux state detected 
by VERITAS. This represents the first detection of a blazar 
flare with VERITAS during bright moonlight observations. 

Within this campaign, no significant variability is detected 
at VHE, nor at lower energies (X-rays and optical). The 
VERITAS light curve is consistent with a constant flux 
(6.3% Crab above 250 GeV, roughly five times the archival 
MAGIC detection) between May 01, 2013 and May 07, 2013; 
additional observations during May 18, 2013 indicate (at a 
2.4cr level) that the high-flux state may have ended at some 
point after the last VERITAS detection on May 07. 

The quasi-simultaneous SED is fitted by a standard 
one-zone SSC model. Even though the nondetection by 
Termi-LAT did not enable a full study of the parameter space, 
the measurements are fully consistent with particle accelera¬ 
tion by relativistic diffusive shocks and simple synchrotron 
and inverse Compton cooling, resulting in a power law with 
injection index around 2.2 and a spectral break of 1 . 0 . 

In recent years the number of VHE blazars has significantly 
increased thanks to the current generation of Cherenkov tele¬ 
scope arrays such as VERITAS, MAGIC and H.E.S.S.. The 
broadband emission and the rapid variability of blazars 
require prompt simultaneous multiwavelength campaigns to 
fully understand the physics of their emitting region. The new 
VERITAS observing strategy under bright moonlight will be 
particularly useful for blazar science, significantly increasing 
blazar monitoring capabilities at VHE. 

Combining data sets taken under various observing con¬ 
ditions is not a new aspect of VHE gamma-ray data analy¬ 
sis. In the case of VERITAS, lookup tables for the analysis 
are generated from simulations which already assume a wide 
range of different zenith angles, azimuth angles (and hence 
geomagnetic field strengths), sky brightness and local atmo¬ 
spheric conditions. The introduction of a new observing mode 
for VERITAS is therefore not expected to adversely affect the 
scientific interpretation of the results. Exploring new modes 
is certainly worthwhile, as the results of this paper demon¬ 
strate. Bright moonlight observing time is particularly use¬ 
ful for variable VHE sources, and can also be used to ob¬ 
tain deeper exposures on steady sources, especially for targets 
with hard VHE spectra which can be detected despite an in¬ 
creased energy threshold. The merging of data taken under 
multiple observing modes and detector configurations will be 
an important analysis challenge for future large arrays, such 
as the Cherenkov Telescope Array (CTA), particularly with 
regard to deep VHE exposures and extended surveys. 
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